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Abstract. We present the calculation of the impact factor for the photon to quark, antiquark and gluon transition within 
Balitsky's shock-wave formalism. We also rederive the impact factor for photon to quark and antiquark transition. These 
results provide the necessary building blocks for further phenomenological studies of inclusive diffractive deep inelastic 
scattering as well as for two and three jets diffractive production which go beyond approximations discussed in the literature. 
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INTRODUCTION 

One of the major achievements of HERA was the experimental evidence that among the whole set of y*p —> X deep 
inelastic scattering events, almost 10% are diffractive (DDIS), of the form y' p —> XY with a rapidity gap between the 
proton remnants Y and the hadrons X coming from the fragmentation region of the initial virtual photon [1,2]. Diffrac¬ 
tion can be theoretically described according to several approaches, important for phenomenological applications. The 
first approach involves a resolved Pomeron contribution (with a parton distribution function inside the Pomeron), while 
the second one relies on a direct Pomeron contribution involving the coupling of a Pomeron with the diffractive state. 
The diffractive states can be modelled in perturbation theory by a qq pair (for moderate M 2 , where M is the invariant 
mass of the diffractively produced state X) or by higher Fock states as a qqg state for larger values of M 2 . Based on 
such a model, with a two-gluon exchange picture for the Pomeron, a good description of HERA data for diffraction 
could be achieved [3], One of the important features of this approach is that the qq component with a longitudinally 
polarized photon plays a crucial role in the region of small diffractive mass M, although it is a twist-4 contribution. 
In the direct components considered there, the qqg diffractive state has been studied in two particular limits. The first 
one, valid for very large Q 2 , corresponds to a collinear approximation in which the transverse momentum of the gluon 
is assumed to be much smaller than the transverse momentum of the emitter [4]. The second one [5, 6], valid for very 
large M 2 , is based on the assumption of a strong ordering of longitudinal momenta, encountered in BFKL equation [7]. 
Both these approaches were combined in order to describe HERA data for DDIS [8]. 

Based on these very successful developments led at HERA in order to understand the QCD dynamics with diffractive 
events, it would be appropriate to look for similar hard diffractive events at LHC. The idea there is to adapt the concept 
of photoproduction of diffractive jets, which was performed at HERA [9, 10], now with a flux of quasi-real photons 
in ultraperipheral collisions (UPC) [11], relying on the notion of equivalent photon approximation. In both cases, the 
hard scale is provided by the invariant mass of the tagged jets. 

We here report on our computation [12] of the y* —> qqg impact factor at tree level with an arbitrary number 
of r-channel gluons described within the Wilson line formalism, also called QCD shockwave approach [13], As an 
aside, we rederive the y* —> qq impact factor. In particular, the y* —> qqg transition is computed without any soft or 
collinear approximation for the emitted gluon, in contrast with the above mentioned calculations. These results provide 
necessary generalization of building blocks for inclusive DDIS as well as for two- and three-jet diffractive production. 
Since the results we derived can account for an arbitrary number of r-channel gluons, this could allow to include higher 
twist effects which are suspected to be rather important in DDIS for Q 2 < 5 GeV 2 [14], 
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FIGURE 1. Left: diagram for y —> qq transition. Right: the 4 diagrams for y —> qqg transition. 


FORMALISM 


As stated before, we use Balitsky’s shockwave formalism. Its application shows that this method is very powerful 
in determining evolution equations and impact factors at next-to-leading order for inclusive processes [15], at semi- 
inclusive level for /?,-broadening in pA collisions [16] or in the evaluation of the triple Pomeron vertex beyond the 
planar limit [17], when compared with usual methods based on summation of contributions of individual Feynman 
diagrams computed in momentum space. It is an effective way of estimating the effect of multigluon exchange. Its 
formulation in coordinate space makes it natural in view of describing saturation [18]. One introduces Wilson lines as 


Ui = U ?i =U ( zt , rj) = Pexp 



(1) 


The operator /? J( is the external shock-wave field built from slow gluons whose momenta are limited by the longitudinal 
cut-off defined by the rapidity rj 

bq = [ -^e- i P*b-(p) 0 (e r '-\p + \). (2) 

We use the light cone gauge st ■ 112 = 0, with szZ being the sum of the external field b and the quantum field A 

g/P=A^+b^ } b ^ 1 (z) =b~(z + ,z)n .2 = 8 (z + )B(z)ri 2 , (3) 

where B(z) is a profile function. The dipole operator U 12 = ^-tr [ U 1 U,) — 1 will be used extensively. 


IMPACT FACTOR FOR 7 -A qcj TRANSITION 

For qq production one can write, after projection on the color singlet state and subtraction of the non-interacting term 


Mq = N c 


Denoting Z\i = \/x q Xqzf 0 , we get for a longitudinal photon 


F {Pc l ,Pq,k 1 Zl,Z2) a Elu = Q{Pq)Q{plj) 

and for a transverse photon 
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IMPACT FACTOR FOR 7 -rf qqg TRANSITION 

For qqg production, projecting on the color singlet state and subtracting the non-interacting term again, one can write 


M a = N 2 Jdzid,Z 2 dziFi {p q ,Pq,Pg,Z0,Zl,Z2,Z3) a ^ (U 32 +U 13 -U 12 +U 32 U 13 ) 
f N 2 — 1 

+ N c l dzidZ 2 F 2 (p q ,Pq,Pg,Z 0 ,Zl,Z 2 ) a (7) 


2 N c 

The first and the second line of this equation correspond respectively to the two last diagrams of the first line and to 
the second line of diagrams of Fig. 1. For a longitudinally polarized photon, they read 
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while for a transversally polarized photon, we have 
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We denote F 2 (p q ,p q ,Pg,Z 0 ,Zl,Z 2 ) a =F 2 (p q ,Pq,Pg,Z0,Zl,Z2) a +JdZ3Fl (p q ,Pq,Pg,Z0,Zl,Z2,Z3) a ■ 


2- AND 3-GLUON APPROXIMATION 


Let us notice that the dipole operator U ; / is of order g 2 . Hence for only two or three exchanged gluons one can neglect 
the quadrupole term in the amplitude M a and get 


Jdzidz 2 Vn ( N 2 - l) F 2 {zi,Z 2 ) a + JdZ3 {N 2 F\ (zi,Z3,Z2)“ +N 2 F\ (z3,Z2,Zl)“ -Fi (zi,Z2,Z3)“} 


M a — 2 I dzidz 2 U12 


( 12 ) 


For p q = pg = p q = 0, those integrals can be performed analytically. Otherwise they can be expressed as a simple 
convergent integral over [0,1] that can be performed numerically for any future phenomenological study. 



CONCLUSION 


The measurement of dijet production in DDIS was recently performed [19], and a precise comparison of dijet 
versus triple-jet production, which has not been performed yet at HERA [20], would be very useful to get a deeper 
understanding of the QCD mechanism underlying diffraction. Recent investigations of the azimuthal distribution of 
dijets in diffractive photoproduction performed by ZEUS [21] show sign of a possible need for a 2-gluon exchange 
model, which is part of the shock-wave mechanism. Our calculation could be used for phenomenological studies of 
those experimental results. A similar and very complementary study could be performed at LHC with UPC events. 
One should note that getting a full quantitative first principle analysis of this would require an evaluation of virtual 
corrections to the y* —> qq impact factor, which are presently under study [22], 

Diffractive open charm production was measured at HERA [23] and studied in the large M limit based on the direct 
coupling between a Pomeron and a qq or a qqg state, with massive quarks [6], Such a program could also be performed 
at LHC, again based on UPCs and on the extension of the above mentioned impact factors to the case of a massive 
quark. This could be further extended to J/ x ¥ production, which are copiously produced at LHC. 
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